Simulation of adsorbate-induced faceting on curved surfaces. 
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A simple solid-on-solid model, proposed earlier to describe over layer- induced faceting of bcc(lll) 
surface, is applied to faceting of curved surfaces covered by an adsorbate monolayer. Surfaces 
studied in this paper are formed by a part of sphere around the [111] pole. Results of Monte Carlo 
simulations show that the morphology of a faceted surface depends on the annealing temperature. At 
an initial stage the surface around the [111] pole consists of 3-sided pyramids and step- like facets, 
then step-like facets dominate and their number decreases with temperature, finally a single big 
pyramid is formed. It is shown that there is a reversible phase transition at which a faceted surface 
transforms to an almost spherical one. It is found that the temperature of this phase transition is 
an increasing function of the surface curvature. Simulation results show that measurements of high 
temperature properties performed directly and after fast cooling down to a low temperature lead to 
different results. 

PACS numbers: 68.35.Rh, 68.43.De, 64.60. Cn, 68.35.Bs, 68.60.Dv 



I. INTRODUCTION 

It has been recently demonstrated that ultrathin metal 
films induce faceting of bcc(lll) surfaces. Atomically 
rough W(lll) and Mo(lll) surfaces, when covered by 
a single physical monolayer of certain metals (Pd, Rh, 
Ir, Pt, Au) and annealed to T > 700K, undergo massive 
reconstruction from a planar morphology to a microscop- 
ically faceted surfaceiAiMiSii. In most of investigated 
bcc(lll) surfaces, the faceted morphology comprises 3- 
sided pyramids with {211} facets and facet sizes ranging 
from ~ 3 to 100 nm. It has been shown that the metal 
that induces the faceting acts like a surfactant and re- 
mains on the surface during the faceting transformation 7 . 

Suggestions that faceting in these systems is ther- 
modynamically favorable have been confirmed by first- 
principles calculations 7 *^ performed for (111), (211), 
and (110) surfaces of Mo and W. Results show lower- 
ing of the surface energy as surfaces are being covered by 
metal overlayer (e.g. Pd, Pt, Au). Moreover, at coverage 
of one physical monolayer, the bcc(lll) surface becomes 
unstable against faceting, i.e. the system gains energy by 
transition from (111) to {211} orientation. 

Interesting results have been obtained for faceting 
on curved surfaces of W and Mo in FIM (field ion 
microscopy )Afliiiiii and FEM (field emission microscopy) 
experimentsiiiAiii. Initially a spherical surface com- 
prises macroscopic crystal facets. After depositing a 
metal film (e.g. Pd, Pt) or oxygen layer and anneal- 
ing at elevated temperatures, an increase of the macro- 
scopic {110} and {211} crystal facets is observed. An un- 
expected change in the faceted morphology occurs near 
the (111) poles, where step-like {211} microfacets are 
formecUfliii. The number of step-like facets decreases 
with the annealing temperature, and in some cases, a sin- 
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gle {211} pyramid around the [111] pole is observed!!. 
Moreover, occurrence of new facets {123}, {178} in Pt 
on W and Pd on Mo has been reported in papers on 

In theoretical studies of complicated surface problems 
(e.g. roughening transition, surface reconstruction, sur- 
face growth, surface phase transitions) simple models 
like lattice gas models or solid-on-solid (SOS) models are 
applied 16 i 17 i 18 i 19 i 20 i 21 i 22 i 28 i 24 i 25 i 26 . In our earlier pape* 2 ! 
we have proposed a SOS model to study the adsorbate- 
induced faceting of the bcc(lll) crystal surface at con- 
stant coverage. Monte Carlo simulation results show for- 
mation of pyramidal facets in accordance with experi- 
mental observations. Moreover, the model describes a 
reversible phase transition from a faceted surface to a 
disordered (111) surface. Such a phase transition has 
been found in the Pd/Mo(lll) system^. In this paper 
we study faceting of a curved bcc crystal surface by us- 
ing the SOS model. We focus on changes in surface mor- 
phology in the vicinity of the [111] pole, where step-like 
{211} facets were observed in FIM experiments. The pa- 
per is organized as follows. In Sec. |H] the SOS model for 
a curved surface of a bcc crystal is described. Results of 
Monte Carlo simulations, including changes in the sur- 
face morphology with the annealing temperature, and a 
reversible phase transition, are presented in Sec. IIIII 



II. THE SOS MODEL 

To study an adsorbate-induced faceting on a spheri- 
cal surface we use a solid-on-solid model of Ref. 27 that 
describes faceting of bcc(lll) surface. Let us recall the 
main assumptions. Atoms along the closed packed direc- 
tion [111] form columns which are represented in the SOS 
model by their positions on the (111) plane, i = (i x ,i y ), 
and heights, hi. Column positions form a triangular 
lattice (see Fig. [Q which is obtained by projection of 
the bcc crystal lattice on the (111) plane. Each column 
comprises substrate atoms and one adsorbate atom at 
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FIG. 1: Schematic top view of the bcc(lll) surface. Atoms 
from three successive geometrical layers represent positions of 
columns in the SOS model. The Z axis is normal to the (111) 
plane and a x = ay/2/2, a y = 



J 



the top. This implies a constant coverage of 1 physical 
monolayer. There is also a restriction imposed on column 
heights: the nearest neighbour column heights can differ 
only by ±1, ±2. A surface atom represented in the model 
by (i x ,iy,hi) has a position (i x a x ,i y a y , hi<2 z ) in the bcc 
crystal, where a z is the distance between neighbouring 
(111) layers (a z = a-\/3/6) and a is the lattice constant. 
The model Hamiltonian represents surface formation en- 
ergy which depends on column heights hi in the following 
form 



H = \ E | E ~ h n I "!) + Ki5 (\hi - h h |-2)] [2 J 2 S (\h t - h h |) + (2 J 2 + K 2 ) 5 - h h \ - 3)] 

+ J 2 J2{ S (\hi-h h \-2) + S(\h i -h j3 \-4)}\+NJ , (1 



where the sums over ji, j'2, and js represent the sums 
over first, second, and third neighbours of the column 
at site z, respectively. Model parameters Jo, Ji, J2, K\, 
and K2 can be expressed by interaction energies between 
substrate and adsorbate atoms (for details see 27 ). 



A. Initial conditions for spherical surfaces 

We assume that the surface has a initially spherical 
shape determined by the radius R and the angular ra- 
dius (see Fig. The integer values of column heights, 
hi, that minimize a deviation from the perfect spheri- 
cal shape, are chosen. Positions of columns on the (111) 
plane form a circle of radius R x = RsinQ. We choose 
fixed boundary conditions, i.e., outside the circle the 
columns form the (111) face and their heights are frozen. 
Such boundary conditions make calculations much easier 
because the number of interaction constants in Hamil- 
tonian Q can be reduced to one, as will be shown in 
what follows. It is worth noting that we can also study a 
fiat (111) surface by an appropriate choice of initial col- 
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umn heights in the circle. In this way we can compare 
adsorbate-induced faceting of flat and curved surfaces 30 
using the same conditions, i.e. the annealing tempera- 
ture, annealing Monte Carlo time, the number of surface 
atoms, and the boundary conditions. The restriction im- 
posed on column heights of the nearest neighbours (see 
Sec. [0} limits the angular radius to < 20°. An exam- 
ple of a spherical surface with = 20° is shown in Fig. |3[ 
where a (111) facet and {211} facets occur together with 
terraces. There are also smaller facets with higher Miller 
indices. It is worth noting that in all surface images 
shown in this paper, the color of a surface atom repre- 
sents the number of its nearest neighbours with smaller 
heights. Such representation allows an easy identification 
of different facets, e.g. {110} facets are represented by 
value 3 and row, trough of {211} facets by 5,1, respec- 
tively. Investigation of faceting of curved tungsten sur- 
face induced by palladium 10 has revealed that step-like 
microfacets are formed in (111) crystal zones. Hence, the 
present model will be used to study faceting around the 
single [111] pole on a spherical surface. 
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FIG. 2: Schematic (a) side view (b) top view of a spherical 
surface. 




B. Energy of spherical surface 

In Ref™ it has been shown that energies per column 
of ideal faces (111), (211), and (110) are the same for the 
interactions limited to the nearest neighbours. It turns 
out that the same is true for the energy of a spherical sur- 
face. We calculated the number of bonds between nearest 
neighbours: b\ with |A/i| = 1 and 62 with |A/i| = 2 for 
hundreds of angles from the interval [0°,20°] and for 
constant number of columns TV. For most of the angles, 
b\ and 62 are exactly the same as in the case of the ideal 
(111) surface (0 = 0°). For a few angles the values of b\ 
and 62 are a little changed e.g., |A&i|/&i = 3.8 x 10 -4 
for R x = 93a x (N = 23497) and |A&i|/&i = 5.6 x 1(T 5 
for R x = 370a x (N — 327367), and such angles will be 
omitted. On the other hand, permitted changes of a con- 
figuration conserve b\ and 62. Thus the energy of the 
nearest neighbours interactions is conserved in this SOS 
model and it can be neglected (or treated as the reference 
energy). For this reason we can limit the number of the 
model parameters to two (J2, K2). Similarly as in2i, we 
can eliminate one more parameter by choosing J2 as the 
unit of energy. We will work with dimensionless quan- 
tities : energy H = H/J2, temperature T = fe^T/^, 
and interaction constant K = K2/J2 (in what follows, 
the tilde will be omitted). It has been shown that the 
energy of the (211) face is minimal when — 2 < K < 0, 
whereas the (111) surface is stable for K > 0. There- 
fore simulations of a flat SOS model 27 with negative K 
(e.g. K = —1.25) show faceting of a bcc(lll) surface. 
In the present model, energy of a spherical surface E sp 
decreases with for negative K but is still greater than 
the energy of a (211) face, e.g. the energy per column 
for K = -1.25 are: E U i = 7, E sp (6 = 20°) = 6.26 
and £"211 = 5.75. Thus we expect that a spherical sur- 
face may undergo faceting after annealing at an elevated 
temperature. In this paper we performed all calculations 
for K = —1.25 for two reasons. First, the same value was 
already used in faceting of a bcc(lll) surface 27 , so one 
can compare results for two different models. Second, we 
estimated the value of K for the Pd on W(lll) system 
using surface formation energies from Fig. 3 in Ref. 7 and 
the relations between Ehki and the model parameters 27 . 
As the first-principle calculations 7 were performed both 
with the local-density approximation (LDA) and with the 
generalized gradient approximation (GGA), we obtained 
two different values of K: —1.06 with LDA and —1.25 
with GGA. We finally chose the latter because the GGA 
is considered as more reliable than LDA 7 . 



III. MONTE CARLO SIMULATIONS 



FIG. 3: Top view of a spherical surface. The color of a surface 
atom represents the number of its nearest neighbours (NN) 
with smaller heights. The brightest color is for 6 NN and the 
darkest one for NN. 



Properties of the spherical SOS model with constant 
coverage are investigated by Monte Carlo (MC) simu- 
lation in canonical ensemble using standard Metropo- 
lis algorithmS2i2S. A new configuration is generated by 
choosing two lattice sites i and j and changing heights of 
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columns at these sites: (hi, hj) — > (hi — 3, hj + 3). This 
is equivalent to moving a substrate atom from a site i to 
a site j. The new configuration is accepted with proba- 
bility p = min(l, exp(— SE/kT), where SE is the energy 
change of the system generated by transition to a new 
configuration. In the simulations we measure the mean 
energy (H), the heat capacity 

C v = ((H 2 ) - (H) 2 )/k B T 2 , (2) 

and the mean-square width of the surface, 




The latter quantity depends also on the curvature radius, 
therefore we introduce the mean-square radial deviation 




where rj is distance of a surface atom j form the center 
of the sphere. 

A. Faceting of curved surfaces 

In this section we investigate changes of a surface cov- 
ered by a physical monolayer of adsorbate during the 
warming up and cooling down processes. Simulations 
start at a temperature To, with the spherical surface de- 
termined by the radius R and the angle (see section 
III A|) . The warming up process is characterized by grad- 
ual elevation of temperature, Ti = T^_i + AT. At each 
temperature the system is annealed for time r, measured 
in Monte Carlo steps per site (MCS). After reaching the 
maximal temperature, the system starts to cool down by 
changing AT — > — AT. In order to measure the tempera- 
ture dependence of interesting quantities we choose 1000 
configurations in the final stage at each temperature and 
we repeat the warming up - cooling down cycle for 20 
samples. In what follows we will explain why averaging 
over samples is important in this problem. 

Results of simulations carried out for K = —1.25, 
R x 93a x , = 20°, AT 0.03 and r = 2 x 10 5 are 
presented in Figs. 0| - The temperature dependence 
of Sr 2 in the warming up - cooling down cycle clearly 
shows (see Fig. 0} that this process is irreversible. For 
T > 0.3 the mean-square radial deviation increases as 
the temperature is elevated, which indicates that {211} 
facets are formed on initially spherical surface. First we 
observe two types of facets on the surfaces (see Fig. [5f a) 
and Fig. EJ: 3-sided pyramids and pits located near the 
pole, and elongated facets forming steps between the 
{211} faces. The steps are built of two alternating facets, 
e.g. (211) and (121). They have characteristic parallel 
long edges, oriented along (311) directions, that can be 
seen in FIM experiment sifiiAi as bright lines. Note that 



5 




0.0 0.5 1.0 1.5 2.0 
T 

FIG. 4: Temperature dependence of 5r 2 during the warming 
up (red circles) and the cooling down (blue triangles) of the 
system. 

there are also concave edges (dark lines on Fig. EJ which 
are not present in FIM images. Further increase of the 
temperature makes the pyramids to disappear and re- 
duces the number of steps (see for example Fig. EJb)). 
In other words, a single pyramid with multiple edges is 
observed. As temperature is increased, even more, the 
number of these edges is reduced (see Fig.^fc)). At high 
enough temperature, it is possible to observe a single 3- 
sided pyramid on the surface (see Fig. EJd)); sometimes 
this pyramid has a double edge. Above Td ~ 1.31, a 
rapid decline of Sr 2 is observed. This indicates a de- 
faceting transition at which the pyramid disappears (see 
e.g. Fig.[5fa)). Although 5r 2 is very small above T^, the 
surface is not spherical. We observe also that the surface 
energy rapidly increases at T^ and the heat capacity at- 
tains its maximum. A similar defaceting transition was 
obtained in simulations of faceting of a bcc(lll) surfaced 
and observed in LEED experiments^ for Pd on Mo(lll). 
The reverse process of slow cooling down the system from 
the highest temperature, leads again to faceting at the 
temperature T^. Hence this is a reversible phase transi- 
tion. 



B. Measurements after fast cooling 

In FIM experiment 11 the number of parallel edges on a 
faceted curved surface was measured as a function of the 
annealed temperature. Because FIM images are taken 
after a rapid cooling to a low temperature, we adopted 
a similar procedure to measure the average number of 
parallel edges, n e . Configurations chosen at a tempera- 
ture T were quickly cooled down (the system was held 
for 20 MCS at each T^) to the temperature T m = 0.2 at 
which calculations where performed. This procedure re- 



FIG. 5: Top view of the surface after annealing to: (a) T = 0.41, (b) T = 0.68, (c) T = 0.86 and (d) T = 1.10. 



duces the thermal disorder and makes it possible to use 
a simple algorithm for edges searching. It is interesting 
to note that the number of parallel edges decreases with 
the temperature down to the defaceting temperature and 
then suddenly increases (see Fig. (TJ. Similar behavior 
was observed in the FIM experiment 11 . The reason why 
fast cooling can reproduce a pyramid with multiple edges 
lays in morphology of disordered surface. When the tem- 
perature is elevated over X^, the pyramid is initially dam- 
aged near the apex and along the edges. In these regions 
the surface is not spherical but has a disordered hill-and- 
valley structure, i.e. there are many small {211} facets 
which form irregular, short {211} steps. As the tem- 
perature is elevated, the disordered regions broaden (see 
Fig. I3a)). On the other hand, fast cooling causes the 
system to return to a pyramid with multiple edges (see 
Fig. E^b)) and multiplicity depends on the width of the 
disordered region and the cooling rate. This fact can be 
used to estimate Td from experimental results, provided 
that the temperature obtained is below the desorption 



temperature. We have also checked how this effect de- 
pends on the cooling rate. If the cooling is 10 times faster, 
then this effect does not occur, i.e. the disordered sur- 
faces remains frozen. In contrast to this, cooling of the 
disordered phase will result in a big 3-sided pyramid on 
the surface for cooling rates over hundred times slower. 

The quantities like n e , £r 2 , 5h 2 have different temper- 
ature behavior for a single sample and when they are 
averaged over many samples. The single sample tem- 
perature dependence of these quantities have step-like 
character (e.g. see Fig. EJ). We think that this effect is 
connected with large sizes of step-like {211} facets on a 
curved surface. A massive transport of atoms is needed 
to reduce one convex edge by merging two neighbouring 
{211} steps. Of course, there should be an energy bar- 
rier for such process and its height must depend on the 
size of {211} facets. Therefore the width of the plateau 
occurring in n e (T) in Fig. increases with the annealing 
temperature. On the other hand, n e (T) averaged over 
many sample becomes a smooth function of temperature 
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FIG. 7: The average number of parallel edges versus temper- 
ature calculated after fast cooling to T = 0.2. 



because the growth of facets proceeds in individual way 
for each sample due to the fact that the distances be- 
tween edges are not constant at a given T (see Fig. 0), 
but they have visible dispersion. 



C. Defaceting of curved surfaces 

It has been mentioned in Sec IIII Al that defaceting 
temperature, X^, of spherical surface with (0 = 20°) is 
greater than in the case of a flat surface. To investi- 
gate the dependence of Td on the surface curvature we 




FIG. 8: Snapshot of the surface (a) at T = 1.4, (b) after fast 
cooling from T = 1.4 to T = 0.2. 

performed a series of simulations for several angles 0, 
keeping the number of columns TV (or R x ) fixed. This 
means that the radius R = R x / sin of a spherical surface 
was varied. Simulation started from a high temperature 
and as the system was cooled down with AT = 0.01, 
we measured the surface energy, the heat capacity, the 
mean-square width and the mean-square radial devia- 
tion near the reversible phase transition. The system 
was held at each temperature for 10 6 MCS. We see in 
Fig HnT a) that the surface energy of a fiat surface rapidly 
decreases at the defaceting temperature T^. This is the 
first-order phase transition and the energy difference is 
equal to the latent heat. Behaviors of E and Sh 2 indi- 
cate that this phase transition differs from the roughening 
phase transitionSii24i25i2& which belongs to the universal- 
ity class of Kosterlitz-Thouless transition. 
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FIG. 9: The number of parallel edges for a single sample (cir- 
cles) and averaged over 20 samples (the thick line represents 
data depicted in Fig[7|. 



It is clearly shown in Fig[IT3|that the temperature of 
defaceting transition increases as the curvature of the sur- 
face increases. For a surface with = 20° the change of 
Td is about 10% with respect to the flat case. We also see 
that the change of the surface energy at the phase tran- 
sition decreases with 0. The same behavior is observed 
in the dependence of Sh 2 on T and (see Fig[HJb)). 

We think that the increase of Td with is related to 
the entropy reduction in the disordered phase. The sur- 
face free energy (SFE), F = E - TS, of the faceted 
phase Ff at the transition temperature Td is equal to 
SFE of the disordered phase Fd- Therefore, the temper- 
ature of phase transition can be obtained from Td(0) = 
AE{0,T d )/AS{0,T d ), where AX = X d - X f for X = 
E, S. As is seen in Fig. ITUT a) the change of internal en- 
ergy AE decreases with at Td. Thus the change of 
entropy AS should also be a decreasing function of 0. 
The decrease of entropy in disordered phase with can 
be justified by the fact that one can observe an increase 
of {211} areas and a contraction of disordered regions on 
the surface at Td as curvature (or 0) increases. 



IV. DISCUSSION 

It is shown that a simple solid-on-solid model can be 
applied to study adsorbate-induced faceting of the [111] 
pole region of a spherical surface. Just as in the case 
of faceting of a bcc(lll) surface, it is possible to specify 
interaction constants that describes a situation when an 
initially prepared spherical surface covered by a physi- 
cal monolayer of adsorbate has a surface formation en- 
ergy greater than {211} surfaces. Thus, the faceting of a 
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FIG. 10: Dependence of (a) surface energy, (b) square mean 
width on temperature near faceting - defaceting transition for 
several angles 0. 



curved surface is thermodynamically favorable. Results 
of Monte Carlo simulations show that the morphology of 
a faceted surface depends on the annealing temperature. 
One can distinguish three stages: (1) Formation of 3- 
sided pyramids around the [111] pole and step-like {211} 
facets between the faces of {211} orientations. The latter 
form long edges observed in FIM experimentsiSiii . (2) 
Disappearing of pyramids around the [111] pole and dom- 
ination of step-like facets whose number decreases with 
the temperature. (3) Formation of a single large pyra- 
mid on the surface. Note that the last two cases were 
observed in FIM experiments. The first case, a mixture 
of two kinds of facets, has not been experimentally de- 
tected yet. Using our model, we can investigate in detail 
the structure of a faceted spherical surface. In particular, 
it is possible to examine properties difficult for experi- 
mental studies: concave edges between elongated facets, 
transition from step-like facets to pyramids and pits, de- 
fects, etc. 

It is shown that there exists a temperature of faceting 
- defaceting transition Td above which disordered regions 
appear along pyramid edges. Although mean-square ra- 
dial deviation 5r 2 rapidly decreases above X^, the surface 
is not spherical. Antczak et al 14 observed in a FEM ex- 
periment of Pt on W tip emitter that faceting diminishes 
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at T > 1600K and the hemispherical form of the emit- 
ter is recovered. As platinum is desorbed from tungsten 
emitter at temperatures above 1900K 14 , this observation 
can indicate that defaceting occurs in Pt on W system. 
It would be interesting to verify experimentally this sug- 
gestion. We found that at temperatures just above X^, 
disordered regions are not " spherical" but they have dis- 
ordered hill- and valley structure with many small {211} 
facets which form irregular, short {211} steps. This fact 
seems to be important for measurements of high- tem- 
perature properties performed after fast cooling down 
to a low temperature, a situation typical e.g. in FIM 
experiments. It is shown that fast cooling may trans- 
form the disordered surface into a multi-step surface (see 
FigJHfb)). This is manifested in a temperature behavior 
of the number of parallel edges n e , which exhibits an 
unexpected increase above the phase transition temper- 
ature Td (see FigEJ). Similar behavior of n e has been 
found for oxygen adsorbed on curved W— . We believe 
that this fact can be used to determine the temperature of 
defaceting, provided that Td is below the desorption tem- 
perature (desorption is not included in the model). For- 
mation of facets after a rapid cooling of disordered surface 
can be explained by a fast massive atomic rearrangement. 
Such effect has been observed in a Pd/Mo(lll) system^ 



in which 3-sided pyramidal facets occur on the bcc(lll) 
surface after annealing to higher temperatures. Above 
Td « 870K the pyramids disappear and surface becomes 
flat. Attempts to freeze a high - temperature flat phase 
ended up with a faceted surface even for 60 K/s of cooling 
rate. 

Another interesting property of faceting found in the 
present MC simulation is the dependence of Td on the 
surface curvature. The temperature of transition Td on 
a fiat surface is smaller than Td on a curved surface. In 
FEM experiment 14 , faceting of a Mo emitter covered by 
Pd has been observed in temperature range 750 — 1040K, 
so at temperatures above Td ~ 870K for the Pd/Mo(lll) 
system^. However, this can not be regarded as the evi- 
dence for higher defaceting temperature in a curved sys- 
tem because FEM observations were made at room tem- 
perature and the cooling effect could affect these FEM 
results. 
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